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Abstract 

We present explicit analytic results for the two- loop top/stop/gluino contributions to the cross 
section for the production of CP-even Higgs bosons via gluon fusion in the MSSM, under the 
approximation of neglecting the Higgs boson mass with respect to the masses of the particles 
circulating in the loops. The results are obtained employing the low-energy theorem for Higgs 
interactions adapted to the case of particle mixing. We discuss the validity of the approximation 
used by computing the first-order correction in an expansion in powers of the Higgs boson mass. We 
find that, for the lightest CP-even Higgs boson, the gluino contribution is very well approximated 
by the result obtained in the limit of vanishing Higgs mass. As a byproduct of our calculation, we 
provide results for the two-loop QCD contributions to the photonic Higgs decay. 



e-mail : 

degrassiOf is .uniroma3 . it 
Pietro . Slavich@cern . ch 



1 Introduction 



One of the main goals of the present experimental program at the Tevatron and at the Large Hadron 
Collider (LHC) is the search for the Higgs boson(s) in order to elucidate the mechanism of electroweak 
symmetry breaking. To support this goal, an accurate theoretical knowledge of the Higgs production 
cross-section, its decay modes, and the important background processes is required (for a general 
review see ref. [1]). 

In the Standard Model (SM) the minimal realization of the Higgs sector predicts a single neutral 
Higgs boson H SM , whose mass can be constrained by electroweak precision data and the direct-search 
limit form LEP to be lighter than roughly 200 GeV. At LHC the main production mechanism for H SM 
is the loop- induced gluon fusion mechanism [2], gg — > H SM , where the coupling of the gluons to the 
Higgs is mediated by loops of colored fermions, primarily the top quark. The knowledge of this process 
in the SM includes the full next-to-leading order (NLO) QCD corrections [3, 4], the next-to-next-to- 
leading order (NNLO) QCD corrections in the limit of infinite top mass [5], soft-gluon resummation 
effects [6], an estimate of the next-to-next-to-next-to- leading order (NNNLO) QCD effects [7] and also 
the first-order electroweak corrections [8, 9]. 

The Minimal Supersymmetric extension of the Standard Model, or MSSM, features a richer Higgs 
spectrum that consists of two neutral CP-even bosons h, H, one neutral CP-odd boson A and two 
charged scalars H^ 1 . As in the SM, the gluon- fusion process is the main production mechanism for 
the neutral Higgs bosons. Furthermore, in the MSSM, the coupling of the gluons to the Higgs bosons 
is mediated not only by colored fermions but also by their supersymmetric partners. Thus, the study 
of the NLO QCD corrections to the cross section for Higgs boson production via gluon fusion in the 
MSSM requires the investigation of a larger variety of diagrams with respect to the SM case. Two-loop 
diagrams involving squarks and gluons were first considered in ref. [10] under the assumption that the 
Higgs boson mass can be neglected w.r.t. the masses of the particles running in the loops. Later this 
approximation has been relaxed and the full dependence on all the relevant masses has been retained 
[11, 12, 13]. Explicit two-loop results for the diagrams that involve quarks, squarks and gluinos, as well 
as for the diagrams that involve the D-term-induced quartic squark couplings, have been presented 
so far only in the limit of vanishing Higgs mass. Indeed, ref. [14] provided analytic formulae for 
the contributions involving top quark and/or stop squarks, valid for vanishing Higgs mass and under 
the additional simplifying limits of zero stop mixing and hierarchical patterns of soft SUSY-breaking 
masses. Ref. [15] presented a computation of the same contributions valid for arbitrary parameters in 
the stop sector. However, the explicit results were too long to be printed, and were made available 
only in the form of a computer code. Recently, the rather daunting calculation of the full two-loop 
QCD Higgs-gluon-gluon amplitude in the MSSM for arbitrary Higgs mass has been completed [16]. 
The calculation relies on a combination of analytical and numerical methods, and explicit results have 
not been made available so far. 

Given the importance of the Higgs-physics search program, it is highly desirable to have the NLO 
radiative corrections to the gluon-fusion cross section expressed in a simple and flexible analytic form 
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that can be easily used to investigate large regions of the MSSM parameter space. To this purpose, we 
present explicit analytic results for the two- loop top/stop/gluino contributions to the cross section for 
CP-even Higgs boson production, under the approximation of neglecting the Higgs mass with respect 
to the masses of the particles circulating into the loops. Our results were obtained by adapting the 
low-energy theorem for Higgs interactions [17] to the case of particle mixing. We confirm the result of 
ref. [15] by a numerical comparison, after accounting for the different renormalization scheme adopted 
in that paper. We then discuss the validity of the approximation of vanishing Higgs boson mass by 
computing the first correction to it, i.e. terms of 0(m 2 l /M 2 ) where M is the generic mass of the 
particles circulating in the loop. We show that, for the lightest Higgs boson h, the gluino contribution 
is very well approximated by the result obtained in the limit of vanishing Higgs mass. As a byproduct 
of our calculation, we provide results for the two-loop QCD contributions to the photonic Higgs decay. 

The paper is organized as follows: in section 2 we summarize general results on the cross section 
for Higgs boson production via gluon fusion. In section 3 we describe our calculation of the two-loop 
contributions to the CP-even Higgs-gluon-gluon form factors. In section 4 we assess the validity of 
the approximation of vanishing Higgs mass. In section 5 we explain how to adapt our results to the 
calculation of the gluonic and photonic decay widths of the Higgs bosons. In section 6 we discuss the 
applicability to other processes of our way to compute the SUSY-QCD corrections. Finally, in the 
appendix we provide the explicit analytical results for the two- loop top/stop/gluino contributions to 
the Higgs-gluon-gluon form factors. 



2 Higgs production via gluon fusion at NLO in the MSSM 

In this section we summarize some general results on Higgs boson production via gluon fusion. The 
hadronic cross section for Higgs boson production at center-of-mass energy y/s can be written as 

a (hi + h 2 ->• <f> + X) = V [ dxidx 2 f a ,hi(xi, p F ) fb,h 2 (x2, Vf) x / dz 5 ( z ^- ) a ab (z) , (1) 

a,b Jo ' Jo V XlX2j 

where 4> = (h,H), = m^/s, fj> F is the factorization scale, f a ,h z ( x , ^f), the parton density of the 
colliding hadron hi for the parton of type a, (a = g, q, q) and a ab the cross section for the partonic 
subprocess ab — > (f> + X at the center-of-mass energy s = xi x 2 s = rri^/z. The latter can be written 
in terms of the leading-order (LO) contribution as 

a ab (z) = a^zG ab (z). (2) 

We consider now the production of the lightest CP-even Higgs boson, h, through gluon fusion. 
The LO term can be written as 



a 



(0) _ "s(^fl) 



T F (-sina-H" + cosa^") 2 , (3) 



128 v 7 ^ 

where G M is the muon decay constant, a s (fi R ) is the strong gauge coupling expressed in the MS 
renormalization scheme at the scale fi R , Tp = 1/2 is a color factor, and a is the mixing angle in 
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the CP-even Higgs sector of the MSSM. Hi (i = 1,2) are the form factors for the coupling of the 
neutral, CP-even component of the Higgs doublet Hi with two gluons, which we decompose in one- 
and two-loop parts as 

Hi = Hf + ^Hf + 0(a 2 ). (4) 

TT 

The coefficient function G a j,(z) in eq. (2) can in turn be decomposed, up to NLO terms, as 



G ab (z) = G^(z) + ^ + 0(a 2 s ), 



TT 



with the LO contribution given only by the gluon-fusion channel: 

: 6(1 - z) S a g S bg . 



The NLO terms include, besides the gg channel, also the one-loop induced processes gq 
qq ->• gh: 



(5) 

(6) 
qh and 



6(1 



TT 



/4 



— sin a H 21 + cos a H\ 



21 



sin a H\ l + cos a Hr, 



it 



+ h.c. 



+ P gg (z) In ( 4") + C A - (1 - z + z 2 ) 2 V x (z) + C A K gg , 



where the LO Altarelli-Parisi splitting functions are 

P„(z) = 2C A 



l„ ( l- 2 ) + Iln(l- 



qg > 



1 



fo(«) + - - 2 + z(l - 2) 



P OT (^) = C F 



1 + (1 



(7) 



(8) 



(9) 



In the equations above, = iV c and = (N 2 — 1)/(2N C ) (N c being the number of colors), 
Po = (11 Ca — 2 Nf)/6 (Nf being the number of active flavors) is the one-loop /3-function of the strong 
coupling in the SM, and 

'ln l (l -z) 



Vi(z) 



l-z 



(10) 



The ^-channel contribution, eq. (7), involves two-loop virtual corrections to gg — ^ h and one- 
loop real corrections from gg — > hg. The former, regularized by the infrared singular part of the real 
emission cross section, are displayed in the first row of eq. (7). The second row contains the non-singular 
contribution from the real gluon emission in the gluon fusion process. The function lZ gg is obtained 
from one- loop diagrams where only quarks or squarks circulate into the loop, and in the limit in which 
the Higgs boson is much lighter than the particles in the loop it goes to 1Z gg — > —11(1 — z)^/(6z). 
Similarly, the functions 7Z qq - and lZ qg in eq. (8) describe the qq —> hg annihilation channel and the 
quark-gluon scattering channel, respectively. They are obtained from one-loop quark and squark 
diagrams, and in the light-Higgs limit they go to lZ q q — > 32 (1 — z) 3 /(27z), 1Z qg — > 2z/3 — (1 — z) 2 /z. 
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The functions TZ. 



Ml 



TZgg, lZ qg are actually completely known. Their expressions can be obtained from 



the results of ref. [18] (see also refs. [19, 20]). 

The one- loop form factors T-L\ and l-Of contain contributions from diagrams involving quarks or 
squarks. The two-loop form factors T-L\ and contain contributions from diagrams involving quarks, 
squarks, gluons and gluinos. Focusing on the contributions involving the third-generation quarks and 
squarks, and exploiting the structure of the Higgs-quark-quark and Higgs-squark-squark couplings, we 
can write to all orders in the strong interactions 



u 2 = 



A* 



m t fi s 2 8 t F t + m\ s 2 p A + X b m b A b s 2 e b F b + 2 m 2 + 2 m| c 2 D b 



X b m b n s 2 e b F b - m 2 z s 2 p D b 



m t A t s 2dt F t +2m 2 t G t 



2 m| si D t 



(11) 
(12) 



In the equations above Xt = l/sin/3 and X b = 1/ cos /3, where tan /3 = v 2 /v\ is defined as the ratio of 
the vacuum expectation values (vev) of the neutral components of the two Higgs doublets. Also, fi is 
the higgsino mass parameter in the MSSM superpotential, A q (for q = t,b) are the soft SUSY-breaking 
Higgs-squark-squark couplings and 6 q are the left-right squark mixing angles (here and thereafter we 
use the notation = sint/?, c v = cos^ for a generic angle if). The functions F q and G q appearing 
in eqs. (11) and (12) denote the contributions controlled by the third-generation Yukawa couplings, 
while D q denotes the contribution controlled by the electroweak, D-term-induced Higgs-squark-squark 
couplings. The latter can be decomposed as 



D, 



hq ^ . ( I'. 



3q 



~ Qq s e w 



q ' 



(13) 



where I^ q denotes the third component of the electroweak isospin of the quark q, Q q is the electric 
charge and 6 W is the Weinberg angle. 

The one-loop functions entering *H\ are 



-iU 



= F, 



G 1 / = 



G 



it 



1 



Ul qi 



Ql\r~ q2 ) 



rn 



2 Go e ( T qi) + _ 2 



qi 



m~ 2 



<3»^ q2 ) + — 2 G\%{Tq) 



rn 



u 



m 



'12 I 



92 



where = 4m^/m^, and the functions Qq 1 and Q\j 2 read 



(14) 
(15) 
(16) 



Gl\r) = r 



Ql%{r) = -2r 



1 + - log 2 



l - 



log 2 



VT 



(17) 
(18) 



The analytic continuations are obtained with the replacement m\ — > m\ + ie . For later convenience, 
we remark that in the limit in which the Higgs boson mass is much smaller than the mass of the 
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particle running in the loop, i.e. r>l, the functions Q\ and Qy 2 behave as 

s ^-l-ik + o ^- <«w-i-£ + o <'- , >- < i9 > 

The discussion above has focused on the production of the lightest CP-even Higgs mass-eigenstate 
h, whose mass is bounded at tree-level by the Z-boson mass and can hardly exceed 130-140 GeV when 
radiative corrections are taken into account (see, e.g., ref. [21]). In this mass range, it can be expected 
that the exact values of the functions and are well approximated by the results computed in 
the limit in which the h boson mass is neglected w.r.t. the masses of the particles running in the loops, 
which we are going to call the vanishing Higgs-mass limit (VHML). In the next section we derive 
explicit analytic results for the two-loop top/stop contributions to T-L\ l and in the VHML, and in 
section 4 we further elaborate on the validity of the approximation. For what concerns the heaviest 
eigenstate H, general formulae for the production cross-section can be obtained straightforwardly 
with the replacements (—sin a — > cos a, cos a — > sin a) in eqs. (3) and (7). However, depending on 
the choice of MSSM parameters, it might not be possible to rely on the assumption that H is much 
lighter than the particles running in the loops. 

3 Explicit two-loop results in the vanishing Higgs-mass limit 

In this section we derive explicit and (relatively) compact formulae, valid for arbitrary parameters in 
the stop sector, for the contributions of two- loop diagrams involving top quarks and/or stop squarks 
to the form factors for the interaction of a CP-even Higgs boson with two gluons in the VHML. 
Besides providing an independent check of the results of ref. [15], our formulae can be easily modified 
to allow for different renormalization schemes for the input parameters. Furthermore we discuss how 
our results for the contributions involving the stop squarks can be adapted to the squarks of other 
flavors. 

3.1 Derivation of the two-loop top/stop contributions 

The starting point of our derivation is the low-energy theorem (LET) for Higgs interactions [17] 
(see also ref. [22]), relating the amplitude M(X,<j>) for a generic particle configuration X plus an 
external Higgs boson (j) of vanishing momentum to the corresponding amplitude without the external 
Higgs boson, M(X). The LET can be stated as follows: the amplitude M(X,<p) can be obtained 
by considering Ai{X) as a field-dependent quantity via the dependence of the relevant parameters 
(masses and mixing angles) on <f>. The first term in the expansion of A4(X) in the Higgs field, evaluated 
at the minimum of the Higgs potential, corresponds to M.(X,4>). Strictly speaking, in case M(X) 
contains infrared (IR) divergent terms the theorem applies to the IR-safe part of the two amplitudes. 
If <p represents a pseudoscalar Higgs boson and X a pair of vector bosons, an additional contribution 
to M (X, cj)) is induced by the axial-current anomaly. This contribution cannot be expressed in terms 
of derivatives of M(X) and must be computed explicitly. 
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To derive the CP-even Higgs-gluon-gluon form factors in the VHML we apply the LET, identifying 
M(X) with the gluon self-energy in the background-field gauge [23]. Then, the top/stop contributions 
to the form factors Hi (i = 1, 2) are given by 

.top/stop _ 2ttv an'(o) 

where v = (vf + v 2 ) 1//2 ~ 246 GeV is the electroweak symmetry breaking parameter, Si (i = 1,2) are 
the CP-even parts of the neutral components of the two MSSM Higgs doublets and II*(g 2 ) denotes 
the top/stop contribution to the transverse part of the adimensional (i.e. divided by q 2 ) self-energy 
of the gluon. In analogy with the discussion in ref. [24], the dependence of the latter on the Higgs 
fields Si can be identified through the field dependence of the top and stop masses and of the stop 
mixing angle. The self-energy depends also upon a fifth field-dependent parameter related to the phase 
difference between the top and stop fields. However, this parameter is relevant only when one consider 
derivatives with respect to the CP-odd fields, thus it can be ignored in our case. As in ref. [24], a 
lengthy but straightforward application of the chain rule allows us to express the functions Ft, Gt , Ft 
and Gt appearing in eqs. (11, 12) and in eq. (13) as combinations of the derivatives of the gluon 
self-energy with respect to the top and stop masses and to the stop mixing angle. In particular, we 
find for the two-loop parts of the functions 

F 2i = _dZ dZ 4c| et dZ ^ 

1 dm 2 dm 2 m? — m~ dcL 

Qit = d^ + d^ + d^^ (22) 

* dm 2 dm 2 dm 2 ' 

ti ti L 

~ 2e = d_Z dZ | As 2 dt dZ ^ 

1 dm 2 dm 2 m~ — m~ 5cL ' 

ti ti ti t 2 2 <>t 

_ az az 

Gt ~ dm[ + dri> (24) 

tl *2 

where, to reduce clutter, we used the shortcut Z = (2/Tp) n 2£, *(0), after decomposing the gluon 
self-energy in one- and two-loop parts as 

Ii{q 2 ) = ^V) + (^Vn 2 V) + 0{0? s ) . (25) 

We computed the contributions to the gluon self-energy from the two-loop diagrams that involve 
top and/or stops with the help of FeynArts [25], using a version of the MSSM model file adapted to the 
background field gauge. After isolating the transverse part of the self-energy with a suitable projector, 
we Taylor-expanded it in powers of the squared external momentum q 2 . The zeroth-order term of 
the expansion vanishes as a consequence of gauge invariance, while the first-order term corresponds 
indeed to n 2£, '(0). We evaluated the two- loop vacuum integrals using the results of ref. [26]. Finally, 
we computed all the derivatives 1 of Z that enter eqs. (21)-(24). 



1 In ref. [26] the two-loop vacuum integrals are expressed in terms of a function $(mf , m|, ra§), whose derivatives can 
be easily obtained using the results of appendix A of ref. [27] . 
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We performed the two-loop computation using dimensional regularization (DREG) and modified 
minimal subtraction (MS). However, it is convenient to express our results in terms of parameters 
renormalized in the DR scheme, which is based on dimensional reduction (DRED) and preserves 
the supersymmetric Ward identities and relations. The conversion of the parameters from the MS 
scheme to the DR scheme was discussed in ref. [28]. In particular, the DR Higgs-quark-quark Yukawa 
couplings differ from their MS counterparts by a finite one-loop shift which, when inserted in the 
one-loop part of a calculation, induces an additional two-loop contribution. On the other hand, the 
couplings of the Higgs bosons to squarks, as far as strong corrections are concerned, are the same in 
both schemes, and they are related by supersymmetry to the corresponding DR Yukawa couplings. 
Specializing to our calculation, only the top contribution to T-L^ is going to induce an additional 
two-loop contribution when the top Yukawa coupling is converted from its MS value to its DR value, 
while the stop contributions to %Y and l-Of can be directly identified as expressed in terms of DR 
parameters. However, as can be seen from eqs. (12, 15, 19), in the VHML the top-quark contribution 
to ti^ goes to a constant, i.e. it does not actually depend on the top Yukawa coupling. Therefore we 
need not introduce any additional contribution 2 to the two-loop results obtained using DREG. 

In the explicit formulae for the derivatives of Z we identify the contributions of diagrams with 
gluons (g), with strong, D-term-induced quartic stop couplings {At), and with gluinos (g). Assuming 
that the parameters in 7i\ e and Ti\ l are expressed in the DR scheme, the non-vanishing contributions 
of the two-loop diagrams with gluons read 



8Z9 
dm 2 



1 



2 to 2 



C F - 



5C A 



(26) 



dm 2 

ti 



2 m? 



3Cf Ca 
4 6 



(27) 



where i = 1,2. Under the same assumption for the renormalization of the input parameters, the 
non- vanishing contributions of the two-loop diagrams that involve strong quartic stop couplings read 



dZ 



it 



dm 2 



Cf 
' 24 



-29 t 



to? + s 2 



W m t 2 



to 



+ 



4 2 

to- vnt 

tl *2 



i 

<s 2 



mf In ^ - m| In 



Q 2 



(28) 



dZ Al 
dm 2 

12 



Cf 
' 24 



4e t m2 r 



ti 



4 

TO- 

t-2 



+ 



4 2 
TOj TO- 
*2 *1 



/ , m~ 



2 ' 
TO- 
tl 



K ln Q2 



(29) 



dZ Ai 
d 4e t 



Cf 
' 24 



(to? 



ti 



TO? ) 2 

12' 



rn 



ti 



TO? 
t2 



2 2 
m- TO- 

tl t2 



In 



rn 



ti 



TO 



f2 



2 

TOj 
tl 



TO 



t2 



Q 2 



TO? 
tl 



2 

TO- 



(30) 



where Q is the renormalization scale at which the DR parameters in the one-loop form factors are 
expressed. Finally, the contributions of the two-loop diagrams that involve gluinos are somewhat 
longer and we report them in the appendix. 



2 Conversely, ref. [16] shows that, if the Higgs-squark-squark coupling is expressed in terms of the MS Yukawa coupling, 
an additional two-loop contribution must be introduced. 
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3.2 Input parameters and renormalization schemes 

To account for the case in which the parameters are expressed in a renormalization scheme different 
from DR, we just have to shift the parameters appearing in the one-loop part of the form factors, after 
taking the limit of zero Higgs mass in the one- loop functions Q}f and Q\j 2 - Since we are focusing on the 
two-loop QCD corrections we need to provide a renormalization prescription only for the top and stop 
masses, for the stop mixing angle, and for the soft SUSY-breaking trilinear coupling A t . Indicating, 
generically, a quantity in the DR scheme as x DR , and the same quantity in a generic scheme R as x R , 
we can write the one-loop relation as x DR = x R + 6x. Then, if the one-loop form factors are evaluated 
in terms of R quantities, the two-loop functions in eqs. (21)-(24) must be replaced by 



? 2i 



11 



F? + 



Gf + 



7T 

6 a, 



7T 



o~m| 
U_ 

4 

ra- 
ti 

c)m~ 
h_ 

4 

m~ 



+ 



<5m~ 

t_2_ 

4 

m~ 

12 
12 



m 



*2 



bm t 5s 2 e t 



m t 



m t 



mt 



+ 



mt 



m 



(31) 



(32) 



G 



21 



Gf + 



+ 



7T 

6 a, 



7T 

6 a, 



(5m? 

M_ 

4 

ti 

u_ 

4 

nix 



+ 



5m2 

t_2_ 

4 

m~ 

12 
12 



m 




(33) 



(34) 



In addition, the two-loop form factor %2 gets a contribution originating from the shift in A t : 



n,2£ 
n 2 



aj2l 
n 2 



m t S 2 t TT 



1 



1 



Sf3 



6 a, 



5A 



m 



rn 



t ■ 



(35) 



A commonly adopted renormalization scheme for the input parameters is the so-called on-shell 
(OS) scheme, in which the top and stop parameters are related to physical quantities. In the OS 
scheme the top and stop masses are defined as the poles of the corresponding propagators, and the 
shifts w.r.t. the DR scheme are 



bra t = ReS t (m t ) , 



5m~ 
ti 



Re IIii (m?) , 



5m~ 

*2 



Ren 22 (m| ) , 



(36) 



where YlJmt) and ITi(m? ) denote the finite parts of the self-energies of top and stops, respectively, 

ti 

each computed at an external momentum equal to the corresponding particle's mass. Explicit formulae 
for the various shifts in eq. (36) can be found, e.g., in eqs. (B.2)-(B.4) of ref. [24]. For the shift in the 
stop mixing angle several OS definitions are possible. We choose 89 t in such a way that it cancels the 
anti-hermitian part of the stop wave- function renormalization (w.f.r.) matrix, leading to [29] 



in 12 (m|)+n 12 (m|; 

69 t = - - 



2 2 

mf — mf 

tl t2 



(37) 
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where n 12 (g 2 ) denotes the finite part of the off-diagonai self-energy of the stops, and is given in 
eq. (B.7) of ref. [24]. Finally, the trilinear coupling A t is related to the other parameters in the 
top/stop sector by 

2m t (A t + fi cot (3) 

s 2 e t = § 2 ' ^) 

m~ — ra- 
ti t2 

therefore, since fj, and tan/3 do not get any 0{a s ) correction, the shift for A t is not an independent 
quantity and it can be expressed in terms of the other shifts: 

5A t = ( Sm l- 5m i + ^_^i) {At + »cotP). (39) 

We have verified that, in the OS scheme, the shifts in eqs. (31)-(35) cancel the explicit dependence of 
'Wy and W?/ on the renormalization scale Q. 

We compared our results with those of the public computer code evalcsusy.f , which is based on 
the results of ref. [15]. The code provides the one- and two-loop parts of the Wilson coefficient for 
the Higgs-gluon-gluon operator in the effective Lagrangian, see eq. (2.1) and (2.5) of ref. [15], using 
an OS renormalization scheme for the parameters in the top/stop sector. After taking into account 
the different renormalization prescription for the stop mixing angle, the different convention for the 
sign of fi in eq. (38), and an overall multiplicative factor in the normalization of the coefficients, we 
find perfect numerical agreement between our results and those of ref. [15]. However, we would like to 
comment on the renormalization prescription for the stop mixing angle adopted in ref. [15] (see also 
ref. [30]). Their counterterm is given by 

se t = , (40) 

mj — mj 

tl t2 

where qo is an arbitrary external momentum (a free input parameter of evalcsusy.f) chosen to be 
of the order of the stop masses. The divergent part of the counterterm for t is indeed compelled 
to have the form of eq. (37) - with the finite part of the self-energy replaced by the divergent part 
- by the requirement that it cancel the poles of the anti-hermitian part of the stop w.f.r. matrix. 
The renormalization prescription given in eq. (40) fulfills this requirement in the case of the QCD 
corrections, because the divergent part of the 0(a s ) contribution to Yl\2{q 2 ) does not depend on q 2 . 
In general, however, this is not the case, unless q 2 , = (m| + m~ 2 )/2. Therefore, we find it preferable 
to stick to the "symmetrical" prescription for 5d t , eq. (37), which can be more naturally applied to 
other loop corrections. 

3.3 Contributions from squarks of other flavors 

The results presented in the previous subsections are valid in the limit in which the Higgs boson mass 
is negligible with respect to the masses of the particles circulating in the loops. Therefore, care must 
be taken in extending the results derived for the top / stop contributions to the contributions of quarks 
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and squarks of other flavors. In the case of the bottom/sbottom contributions the general formulae of 
section 2 hold. However, the VHML can be strictly applied only to the two-loop contributions arising 
from diagrams with sbottoms and gluons, and to those arising from diagrams with quartic sbottom 
couplings. Results valid in the DR scheme can be obtained for the former from eq. (27) and for the 
latter from eqs. (28)-(30), with the trivial substitution t — > b in the squark masses and mixing angle. 
Obviously, the VHML cannot be applied to the contributions arising from two-loop diagrams with 
bottom quarks and gluons, but the exact results for those contributions are available in the literature 
[4, 11, 12]. For what concerns the contributions of two-loop diagrams with bottom, sbottom and 
gluino, the VHML can only be applied under the further approximation that the bottom mass and 
the left-right mixing in the sbottom sector are set to zero (i.e. m& = 0& = 0), effectively killing the 
Yukawa-induced interactions between Higgs bosons and bottom (s)quarks. Since the left-right sbottom 
mixing contains a term proportional to mj, t an P, this is not a good approximation when tan f3 is large 
enough to offset the smallness of m^. 

When the bottom Yukawa coupling is neglected, the only diagrams that give a contribution to 
the form factors for the Higgs-gluon-gluon interaction are those in which the Higgs boson couples to 
the sbottoms through the electroweak, D-term-induced interaction. These diagrams contribute to the 
function Db that appears in eqs. (11, 12) and is further decomposed into two functions and Gb in 
eq. (13). The expressions in eqs. (23) and (24) for the two-loop part of the functions simplify to: 

m _ _ez dz_ ~ 2e az dz 

b dm? dm? ' b ~ dm? + dm? ' 1 j 

ol or °l or 

where, in the absence of left-right mixing, the sbottom mass eigenstates b\ and b<i are identified with 
bL and 6r, respectively. The DR contributions from the two- loop diagrams with gluons and with 
quartic sbottom coupling can again be read off eqs. (27)-(29) after replacing t — > b and setting 8b = 0. 
The contribution of the two-loop diagram with gluino, sbottom and (massless) bottom reads 3 

dZ'a m? , m? ml ( m? + m? 



dm? 6m~ Q 2 (m? — m?) 2 \ m? 




2+ 9 h 




12 (mi-mi) 



(42) 



where i = L,R. If the sbottom masses appearing in the one-loop part of the form factor are taken 
as the physical ones, the functions F b and G 2 must be shifted as in eqs. (33) and (34), neglecting 
the term proportional to 5c2e b - It is useful to remark that, as is clear from eq. (13), there is a partial 
cancellation between the sbottom contributions in Dj, and the corresponding stop contributions in D t . 
Therefore, the sbottom contributions controlled by the electroweak gauge couplings must be taken 
into account even when tan f3 is small and the bottom Yukawa coupling can be neglected. 



3 Note Added: We later found that eq. (42) is incorrect. The term proportional to Cf in the contribution of the 
two- loop diagram with gluino, sbottom and bottom can in fact be obtained by setting 8 b = in eq. (40) of ref. [38]. 
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Finally, for the squarks of the first two generations the approximation of neglecting the Yukawa 
couplings is always satisfactory. The remaining D-term-induced contributions to the form factors can 
be obtained by trivially adapting the results in eqs. (27)-(29) and (42). 

4 On the validity of the vanishing Higgs-mass limit 

In the previous section we presented analytic results for the two-loop form factors Hj e valid in the 
VHML. As already mentioned, we can expect this approximation to be quite good for the lightest 
CP-even Higgs boson h, while for the heaviest Higgs boson H it is probably less accurate. To put this 
expectation on a more solid ground we computed directly the two- loop top/stop contribution to the 
Higgs-gluon-gluon amplitude via a Taylor expansion in the external Higgs momentum up to terms of 
0(m^/M 2 ) - where <fr = h,H, and M denotes generically the masses of the heavy particles in the loop 
(i.e. top, stops and gluino). The validity of the Taylor expansion is restricted to Higgs masses below 
the first threshold that is encountered in the diagrams. This is always the case for h, while for H this 
situation is realized only in specific regions of the parameter space. 

We computed the Higgs-gluon-gluon amplitude following the same strategy employed for the cal- 
culation of the gluon self-energy, see section 3.1. The zeroth-order term in the Taylor expansion 
reproduces the result that we obtained via the LET, while the 0(m^/M 2 ) term in the expansion gives 
the first correction to the VHML. For simplicity we neglected all the (small) D-term-induced elec- 
troweak contributions. We expressed our results in the OS renormalization scheme outlined in section 
3.2. We remark that, when converting to the OS scheme the results obtained originally in the MS 
scheme, we must introduce additional two-loop contributions of 0(m'^ ) /M 2 ), originating from the shifts 
in the parameters that appear in the Oim^/M 2 ) parts of the one-loop form factors - see eq. (19). We 
checked that the additional contributions cancel the explicit renormalization-scale dependence of the 
0{m'^ ) /M 2 ) part of the two-loop form factors. The analytic expressions for the 0{m 2 ^/M 2 ) corrections 
are very long and we do not report them. 

Similarly to ref. [15], we present our results in two representative MSSM scenarios. The first 
scenario is the so-called SPSla' slope [31], in which the soft SUSY-breaking parameters at the GUT 
scale are related as 

m = 0.28 m l / 2 , A = - 1.2 m 1/2 , (43) 

where mo and are universal SUSY-breaking masses for scalars and gauginos, respectively, Aq is a 
universal Higgs-sfermion-sfermion interaction term, and the other relevant parameters are tan f3 = 10 
and jJL < (with our sign convention). We vary the GUT-scale gaugino mass m 1 / 2 between 100 
GeV and 400 GeV, and the other SUSY-breaking parameters as in eq. (43) above. To compute the 
superparticle masses, we evolve the soft SUSY-breaking parameters down to the weak scale using the 
public computer code SoftSusy [32]. In the second scenario the light Higgs boson is "gluophobic" , 
i.e. the top quark contribution to the Higgs-gluon coupling is largely canceled by the contribution of 
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Figure 1: Two-loop form factors and T-Ljj in the SPSla' scenario (upper plots) and in the gluophobic 
scenario (lower plots). The plots on the right include only the superparticle contributions. In each plot, the 
solid lines refer to the result obtained in the limit of vanishing Higgs mass, while the dashed lines include 
also terms of 0(m2/M 2 ). The lines for 7-Ljj in the SPSla' scenario are truncated where thr ~ 2mf. 
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a light stop. The scenario is defined directly in terms of weak-scale parameters, which we choose as 



7T 



tan/3 = 10, 9 t = j, m A = 300GeV, fi = -500 GeV, m~ g = 500 GeV, = 200 GeV, (44) 

while mi 2 is varied between 250 GeV and 750 GeV. In the first scenario the CP-even Higgs masses 
and mixing angle are computed directly by Sof tSusy, while in the second scenario we compute them 
using the two-loop 0(ata s ) results of ref. [24]. In both scenarios we take nit = 172.6 GeV [33]. 

The four plots in figure 1 show our results for the two-loop form factors Tiff and T~Ljf, which we 
define as 

%f = T F (-smaUf + cos aUf) , H% = T F (cos a %f + sin a nf) . (45) 

The two upper plots refer to the SPSla' scenario, while the two lower plots refer to the gluophobic 
scenario. For each scenario, the plot on the left shows the complete two- loop top/stop contributions 
to the form factors, while the plot on the right shows only the contributions of the diagrams that 
include superparticles. Finally, for each Higgs boson <p = h,H , the solid line corresponds to the result 
obtained for % 2 ^ in the VHML, while the dashed line includes also the contribution of the first-order 
term in the expansion in powers of rn^/M 2 . 

It can be seen from figure 1 that, for the lightest Higgs boson h, the corrections of C(m^/M 2 ) 
to the results obtained for Tiff in the VHML are quite small, which should not come as a surprise 
since m/, is always considerably smaller than m<. In the SPSla' scenario, which as m 1( / 2 increases 
is characterized by relatively heavy superparticles, the comparison between the plots on the left and 
right sides shows that the bulk of the corrections is contained in the diagrams with top quarks and 
gluons, for which a complete analytic result (i.e., valid for any value of the top and Higgs masses) is 
available [4, 11, 12]. In the gluophobic scenario there is still a small Higgs-mass dependence in the 
superparticle contribution, due to the presence of a relatively light stop. In summary, it appears to 
be quite safe to approximate the two-loop top and stop contributions to Tiff with the results obtained 
via the LET in section 3. In case a more refined approximation is required, one can implement the 
complete result for the top/gluon contribution by replacing 



8Z9 



dm 2 2 m 2 



(46) 



in the gluonic part of the function G 21 defined in eq. (22). The functions Q^^ R> and Q^^ Ca ^ in the 
OS renormalization scheme are defined in eqs. (2.15) and (3.8) of ref. [12], respectively. 

The situation is quite different for the heaviest Higgs boson H. As we mentioned before, the limit 
of vanishing ran can only be considered if tuh is smaller than the lowest threshold appearing in the 
loops. Indeed, the curves for T-Ljj in the SPSla' scenario are truncated around m X j 2 = 200 GeV, where 
iuh approaches 2m t . On the other hand, it can be seen from figure 1 that, in the regions of the 
parameter space where the limit of vanishing niH can be applied at all, the resulting approximation 
is not exceedingly bad even for m# as large as 300 GeV. It should however be kept in mind that, in 
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the case of large tan /3, the VHML approximation cannot be applied even if m# is relatively small, 
because the H boson has enhanced couplings to the bottom (s)quarks (see the discussion in section 
3.3). In that case a full computation is unavoidable [16]. 



5 Gluonic and photonic Higgs decays 

The results of the previous sections can be directly applied to the NLO computation of the gluonic 
decay widths of the CP-even Higgs bosons in the MSSM. We specialize to the decay width of the h 
boson, but the formulae presented below can be applied also the decay of H using the replacements 
indicated at the end of section 2 (with the same restrictions outlined there and at the end of the 
previous section). 

At NLO in QCD the decay width of the lightest Higgs boson in two gluons reads 



T(h -> gg) = G » a ^f™l 



Tp (- sin a n\ e + cos a U^) 



7T 



1 + -C ) , (47) 



where C = C v i r t + C ggg + C gq q includes contributions from the two-loop virtual corrections and from 
the one-loop real radiation processes h — > ggg, h — > gqq. The contribution of the two-loop virtual 
corrections is straightforwardly obtained from the results of the previous sections: 

r, „ , ( A\ (- a H"f + cos a H'i> u \ 

C „ rt = C A - + ft 1„ (£| j + + h ,.j , (48) 

while the contributions of the real radiation processes are, in the VHML, 

C 999 = ~ Ca [-j- - , C m = - 7 -N f , (49) 

where Nf is the number of light quark species, with the quarks treated as massless particles. 

Considering for illustration the two scenarios described in section 4, we found that the numerical 
value of the last term (within parentheses) in the r.h.s. of eq. (48) is approximately — 2.6 Cf + 2.5 Ca 
in the SPSla' scenario with = 200 GeV, and — 4.5 Cf + 2.8 Ca in the gluophobic scenario with 
mi 2 = 500 GeV. For the heaviest Higgs boson the corresponding values are —2.9 Cf + 1.4 Ca (SPSla') 
and —1.1 Cf + 0.9 Ca (gluophobic). In both scenarios the superparticle contribution to the coefficient 
of Cf is comparable to the corresponding SM contribution, while the superparticle contribution to 
the coefficient of Ca is much smaller than its SM counterpart. 

As a byproduct of our calculation, we can also provide the explicit results for the two-loop QCD 
corrections to the quark/squark contributions to the photonic Higgs decay. The partial width for the 
decay of the lightest CP-even Higgs boson h in two photons can be written as 



r(/i-> 77 ) = 



G> "em m h 

128V2tt 3 



— sin a V\ + cos a V2 



2 



(50) 
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where a em is the electromagnetic coupling and V\ (i = 1, 2) are defined, in analogy to the Hi in eq. (3), 
as the form factors for the coupling of the neutral, CP-even component of the Higgs doublet Hi with 
two photons. At one loop, the form factors V\ and V2 receive contributions from all the electrically 
charged states of the MSSM (see, e.g., the first paper in ref. [4] for the explicit results) and they are 
in general dominated by the contribution of the diagram involving the W boson. However, only the 
contributions involving quarks and squarks receive QCD corrections at two loops. We separate the 
one-loop part of the form factors and the two-loop QCD corrections as 



u 



+ ^Vf + 

7T 



(51) 



where the ellipses stand for three- loop terms of 0(a 2 s ) and for two-loop terms controlled by other 
coupling constants. 

Focusing on the contributions of the third-generation quarks and squarks, V\ and V2 can be 
decomposed exactly as in eqs. (11) and (12), with the following substitutions in the r.h.s. of eqs. (11) 
and (12) 



9' 



G q -> G q , D q -> D, 



Q 2 t N c 



Ql N c 



(52) 



sin (3 ' cos (3 

where Q q is the electric charge of the quark. With our overall normalization we have for the functions 
entering the one-loop parts of the form factors, Vf e , 



F u = F 
q q 



while for the ones entering V, 



2t 



G\ l = G 



u 



u 



? 2t 



G 



2t 



G 



2(: 



C A =0' 



D 



2i 



D 



2li 



C A =0 



(53) 



(54) 



q q c A =o' 

i.e. the functions entering the two-loop parts of the form factors can be obtained by setting Ca = 
in the results presented in section 3 and in the appendix. 



6 Discussion 

The LET allowed us to derive explicit and compact analytical formulae, for the top/stop contributions 
to the form factors for the interaction of a CP-even Higgs boson with two gluons, valid in the limit in 
which the mass of the Higgs boson is neglected w.r.t. the masses of the particles running in the loops. 
By direct inspection of the first correction to the results obtained in the VHML, we have argued that, 
for the lightest MSSM Higgs boson, the VHML results provide a quite good approximation to the full 
result, whereas for the heaviest Higgs boson the approximation is less good, and can be applied only 
in specific regions of the parameter space. For what concerns the sbottom contributions to the form 
factors, the validity of our results is limited to the case of small or moderate tan/3. 

As mentioned in section 3.1, the form factor for the interaction of a CP-odd Higgs boson A with two 
gluons receives an additional contribution from the axial-current anomaly and cannot be computed 
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in the same way as the form factors for the CP-even bosons (indeed, the first derivatives of the gluon 
self-energy w.r.t. the CP-odd parts of the Higgs fields vanish at the minimum of the Higgs potential). 
In the VHML the A-gluon-gluon form factor does not receive any contribution beyond one loop from 
diagrams involving gluons [34], as a consequence of the Adler-Bardeen theorem [35]. However, there is 
a non- vanishing two-loop contribution from the diagrams involving gluinos which requires an explicit 
diagrammatic calculation [36]. 

Finally, the LET can also be applied to multiple Higgs boson production [22]. By taking multiple 
derivatives of the gluon self-energy w.r.t. the CP-even or CP-odd Higgs fields, it is indeed possible to 
compute the SUSY-QCD corrections to the production of any number of CP-even Higgs bosons and of 
any even number of CP-odd Higgs bosons in the VHML. The simplest case would be the calculation 
of the SUSY-QCD corrections to the pair production of both CP-even and CP-odd Higgs bosons. An 
analysis of these processes, restricted to the top-quark contributions, has been presented in rcf. [37]. 
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Appendix 



We provide in this appendix the explicit expressions for the contributions of the two-loop diagrams 
with top, stop and gluino to the derivatives of Z. For every derivative we separate the coefficients of 
the color factors Cf and Ca as 



dxi 



= C F 



dxi 



+ C A - 



dz g n 



dxi 



(Al) 



for Xi = (mf , m-^ , m-^ , c 2 ,^). In terms of the two-loop function 5>(x,y, z) defined, e.g., in the 

V of ref. [27], and usin 
the contributions to the various derivatives 4 read 



appendix A of ref. [27], and using the shortcut A = m~ + + m\ — 2 (m~ + m| m 2 + m 2 ^ m 2 ) , 



dm 2 
ti 



+ 



+ 



1 



6mf A 2 

ti 



(m~ + m 2 ) A 2 + 2 m~ m 4 (A + 10 mj m 2 ) 



6 mt mf A 2 



mfA 2 + 2 m~ (A + 5 m~ mf) (m~ + mf 



m x 



4 2 

m§ log^f \ (mf. -mf -4m, 2 ) A 2 



6 m 4 A 3 m? 

ti 1 



+ 2 m? (-18 m 2 m? A + ((3 mf - m~) A - 30 m~ m 2 ) (m 2 - m~ + )) 



^2 m fi (m fi +m t ) + m t (rrig — m t )) A 



m^mt 

2 m 2 |-9m- m^ A + ((2 mj — 9 m^ ) A — 30 m~ m, 2 m^ ) (mj — m 2 + m^ )^ 



2 

m~ 



6 A 3 



m 



log— |- \ A 2 + 12 m\mf A + (2m~ A + 60 m~ra 2 m~ )(mjf + m~ - mf ) 



mr 



| m 2 m\)[m 2 + m |-m| i/ 



2g2e t 
m g m t 



(m~ + m 2 ) A 2 + m~ m 2 (3 mj + 3 m 2 + 20 ) A + 60 m t m| 



I/22 \ i m f 

— j- [m g + m t - s 2dt m- g m t ) log 

(mj , m, 2 , ) | (m~ + mf + 3 m^ ) A + 20 m~ to 2 



4 2 



g20 t 

m g m t 



A 2 + 2 m~ m 2 A + (3 m~ A + 10 m? m 2 m| ) (m| + mf - m~ ) 



(A2) 



dZ^ /dm~ 2 can be obtained from <9Z 5 /dm- through the replacements £i — > £2 and s 2 e t — > — s 2e t 
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dm- 
ti 



+ 



+ 



+ 



12 A* P'"' A 




2 m? 


3m- g A 2 


m? 
12 A3 log 


m? 

2 

ra- 
ti 


2 s 2 e t m t r 



2 2^ 



2 2 i 
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^1 



+ 



S20 t 
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+ ((m? + m 2 + 2 ) A + 20 m? m 2 (m 2 — m? + m 

(m? - m\ + m| ) A 2 + 2 m? (2 m? (m? - m| ) A 
+ ((m? + 5 m| ) A + 10 m? m 2 ) (m? — m. 2 + 



(A3) 
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ti 
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4 m? A + ( A - 10 m? mf ) (m? - m 2 - m? 



A 2 + 2 (5 m? m? A + ((2 m? - m?) A + 10 m? m 2 m? ) (mf - m 2 + m? )) 



+ 



4 



2 
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